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ABSTRACT 





In this study, the author uses a MATLAB simulation to 





develop and test a generalization of the traditional Random 





Search model which allows both the searcher and target to 
move and to be in different, but overlapping, areas. Also 
the best evasion speed for a randomly moving target against 


a Systematic Search is studied. 
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EXECUTIVE SUMMARY 





In this study, the author uses a MATLAB simulation to 





develop and test a generalization of the traditional Random 


Search model which allows both the searcher and target to 





move and to be in different, but overlapping, areas. Also 
the best evasion speed for a randomly moving target against 


a Systematic Search is studied. 


The new generalized Random Search formula is, 


F(t) = AF pig (t) + U- @) Fry, (0), Where 


a =exp(-Evr) , 
A,U 


JA,ARV ” 
Fryoy(t) = (A, 1A, —exp(-2RVt/ A,)) , 
Fos: (t) =1—exp(-2RV(A,,/A,)t/A,) , and 


V = mean relative speed between searcher and target. 


With regard to the best evasion speed against a 





Systematic searcher, xtensiv simulation suggests that a 





5%-20% of searcher speed is the best speed for the target. 
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I. INTRODUCTION 


A. BACKGROUND 


Random Search is a popular model for area search 


because it is mathematically simple and it provides a 





conservative, lower bound on the probability of detecting a 








stationary target with Systematic Search. However, the 





Random Search model also has significant limitations. In 


particular, it assumes that the searcher and target are 








contained in the same area and that the target is stationary. 


We address these model limitations in this thesis. 





B. RESEARCH QUESTIONS 


This thesis will address the following two questions: 


1. How can the Random Search model be generalized for 





Situations wher th searcher and target areas are not 





coincident and both the searcher and target are moving? 


2. What is the best speed for a randomly moving target 





to evade a searcher conducting a Systematic Search? 


Cc. THESIS ORGANIZATION 


The thesis is organized in the following manner: 








e Chapter reviews the search models of Exhaustive 





and Random Search. 














e Chapter presents a MATLAB simulation of Random 








Search. 





Chapter IV develops the extended Random Search 


model wher searcher and target areas are not 





coincident and both the searcher and target are 


moving. 


Chapter V investigates the optimal speed for a 


target evading a Systematic Search. 


Chapter VI summarizes the results of simulations 








and recommends future studies. 


II. OVERVIEW OF RANDOM SEARCH THEORY! 


A. EXHAUSTIVE SEARCH 


dy, Characteristics of Exhaustive Search 


Exhaustive Search is an idealized search model which 





assumes a stationary target and a ‘perfect search’, meaning 





no search overlap, no search effort placed outside the 


target area, and all of target area is completely covered by 








the search sensor. The searcher has a ‘cookie-cutter sensor’ 


with range R, sometimes called a definite range law sensor. 





Such a sensor always detects a target within a specified 





range R, and never detects targets beyond that range. That 


is: 


; 1, range<R 
P(detection)=P, = é ‘ (1) 
, range>R. 


Because of its optimistic area coverage assumptions, 





Exhaustive Search is generally assumed to provide an upper 


bound on the performance of realistic search. 


2. Exhaustive Search for a Uniformly Distributed 
Target 


Assume that a searcher conducts Exhaustive Search over 





the search area A(mn’) with a cookie-cutter sensor of range 


R(nm), speed V (nm/hour), and sweep rate 2VR (nm/hour). The 





1 The references for this section are the unpublished lecture notes 
of Professor James N. Eagle and Professor Alan R. Washburn’s ‘Search and 
Detection’ book (2002). 





time of detection in the random variable TJ (hour), which is 





uniformly distributed between time 0 and A/(2VR). That is, 





T~U[0, A/QVR) ]. 


The probability density function and cumulative 


distribution function for random variable T are: 


(2VR)/ A, t €[0, A/(2VR)] 

0, otherwise 

0, <0 (2) 
F(t) =4 (2VRt)/ A, t €[0, A/(2VR)] 

1,t> A/(2VR) 


f(b) -| 


The mean time to detection is given by the expected 


value of T: 





E(T)=A/(AVR) . (3) 
f (t) F(t) 
1 
(2VR) /A 
A/ (2VR) A/ (2VR) 
Figure 1. Density and CDF for the Time of Initial 














Detection JT with Exhaustive Search. 


B. CONTINUOUS SEARCH 
1. Model Assumption and Definitions 
e y(t) = Ceteckion: rate at: time ¢ tuniis: olptame)., 








Detection rate has two equivalent interpretations: 
4 


i. P(detection in [t,t + At]) = y(t)At+ o(Ar) 


‘ (4) 

= v(t)At, for small Ar 
Note: o(At) is a function of At which goes to 0 
faster than linearly. That is: 


ia ee 
At>0 At 


0 


ii. E(# detections in [t,t+At]) = [ras (5) 
= y(t)At, for small Ar 


T is the random time of initial detection, and we 


define F,(t)=P(T <n), and G,(t)=P(T >t) =1-E,(t) . (6) 





Detection events in non-overlapping time intervals 











are probabilistically independent. So, 
P{N(t+ At) =0} = P{N(t) =0, N(t+ At) — N(t) =0} 
= P{N(t) =0}- P{N(t+ At) — N(t) = 0} (7) 


where, N(t) be the random number of event occur during (0, f] . 


Derivation of G,(t), F,(t) and E(T) 


Since the events of detection (and non-detection) are 





independent in non-overlapping time intervals, 
G, (t+ At) =G,(t)d—-[y@)At + o(At)]) . (8) 
KH 
A B 
Where: 








A=P (non-detection from time to 0 to time ft) 

















B=P (non-detection from time t to time f+At). 





Re-arranging terms, 





At) — A 
Or CN SOM ace) , Letting At go to 
At At 
d ; ; 
0, po Oe) . The solution to this differential 
equation is: 
G,(t) =exp(-[7(s)ds) . (9) 
Therefore: 
F,(t) =1-G,(t) =1-exp(-[ 7(s)ds) 
d t 
frlt) =—F,() = y(Hexp(-[ 7(s)ds) (10) 
dt 0 
E(T)= J thr @adt = [ot y(t)exp(— I (v(s)ds)dt 
If [_, frat =1 , then 7T is a proper random variable, 
eventually detection is certain, and a slightly simpler 








expression is possible. For any proper non-negative random 


variable T: 


E(T) =|" (1- F,(a)dt =f" exp(-[. (s)ds)dr 


The proof is sketched below in Figure 2. 





(11) 








fo — F(x))dx 
zr 
area above F(z) and below 1 


~ > aAF(2i) 
> tif (ai) Ax 


vi 


2 


[et@ae = E(X) 


Dji+]1 





Figure 2. E(X) for Proper, Non-Negative Random 
Variables X 


Cc. RANDOM SEARCH 


1: Characteristics of Random Search 


Assume a search where: 





Searcher has a cookie-cutter sensor. 


b. Each small segment of the searcher’s track is 





randomly and uniformly distributed over the search 





area. 


No search effort falls outside the search area. 


d. Target position is fixed in search area. 





During any time interval At, the searched area is 





2RVAt . Since this area is uniformly distributed over the 





search area A, the probability of its covering the target 
is: 


2RVAt ; 
id = P(detection during [f, t+ Ar]) . 





Thus, we have a constant detection rate with y=(2RV)/A. 
Therefore: 
F(t) =1—exp(—2RVt/ A) 


f-(@) =(2RV/ A)exp(-2RVt/ A) . (12) 
E(T) = A/QRV) 


2. Dynamic Enhancement of Random Search 


Now assume the target with speed U and the searcher 
with speed V move randomly and independently over area A. 


The searcher sweep width is 2R, and T is the random time 





of initial detection. As before, we wish to compute 


F.@)=P7 st) , and for U<<V , we would expect that 





F,.(t)x1-e7R""4 , And more generally, we will look for a speed 


V>max(V,U) such that F,(t)x1—e° a4 , 





V (0) = JV? + U2 — 2UV cos(6) 


~~. 





Figure 3. Dynamic Enhancement of Random Search 


Referring to Figure 3, V(@) is the relative speed 





between searcher and target when @ is the difference in 


their courses. By the Law of Cosines, V(0) =JV? +U’ —2UV cos(6) . 


We now assume that 9 is uniformly distributed between 0 and 





22, so the average V(@) is: 


V=(" vi) fide 


6=0 
2a 2 2 1 

= I JV? +U? —2UV cos(0)(—_)d0 
oy 20 


(13) 
_1 2 2 
=o LAV +U~ —2UV cos(@)d0 





> max(V, U). 


This equation suggests that in Random Search with dynamic 





enhancement, the searcher speed is effectively “enhanced” by 








the target speed. The author will confirm this by MATLAB 








simulation in Chapter 

















3. Random Search when Searcher and Target Patrol 
Areas are not Identical 


As : Searcher’s Search Area 
Ast : Overlap Area 


At : Target Movement Area 








Figure 4. Tllustration of ‘Search Area is not 
Completely Overlapped with Target Area’ 














fast target 


slow target 


time, ¢ 


Figure 5. F(t) for Fast and Slow Target Assumption 


10 


As shown above in Figure 4, we now consider the 


Situation where the search area is not completely overlapped 











with the target area. Two types of target behavior were 











assumed in order to get a simplified mathematical expression. 





a. Slow Target 


Assumptions: 





- Target initial position is uniformly distributed 





over Ais 





- A target starting inside (outside) A, will remain 











so for th ntire search time. 





In this case F(t) is: 


F,.(t) = P(detection by time t | target start in A, )P(target starts in A, ) + 
P(detection by time t | target starts outside A,,)P(target starts outside in A,,) 
= P(detection by time t | target starts in A, ) P(target starts in A, ) (14) 
Bs 
=( 5 


t 


)(1—exp(—2RVt/ A,)). 





b. Fast Target 


Assumption: 


- Target spends (A,/A) of the search time inside 





A, . 


Ss 


By time ¢t, this target has been available for 


detection for (A,/A)t time units. Therefore, 


F(t) =1—exp(-2RV(A,/A)t/A,) . (15) 


11 





In light of these assumptions, eventual detection 





for the fast target is assured. In contrast, the upper bound 





of F(t) for the slow target is (A,/A). In this thesis, we 
will attempt to generalize the fast and slow target models 


to allow F,(t) to be estimated for intermediate target speeds. 


12 


IIIT. SIMULATION OF RANDOM SEARCH 


A. DESCRIPTION OF RANDOM SEARCH MODEL 


dy, Characteristics of the Searcher 


The searcher is assumed to have a “cookie-cutter” 





sensor of radius R . In addition, the searcher is not 





allowed to search for a target outside the search area. 


Therefore th searcher’s allowable position is limited by 














sensor radius R. For example, if the search region’s X - 
axis length is 50nm and the Y-axis length is also 50nm, 


then the overall search area (A) is 50x50nm’. But the 








searcher’s actual moving area (A') is (50—R)x(50—R)nm’ . 








Th searcher’s initial position is uniformly 





distributed inside of the search region A'. After that, the 








searcher chooses his course randomly, independent of the 


target’s movement. The course change event is determined by 








Poisson process with rate A. It is assumed that the speed 





of the searcher is always faster than that of the target. 








This is allowed because the searcher and target roles can be 
reversed in the simulation, resulting in the same 


probability of detection. 


2. Characteristics of the Target 


The target’s initial position is uniformly distributed 





over the search region A. The logic of the target movement 


is the same as that of the searcher, and the target has its 


own, independent Poisson process with course change rate 4. 
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B. COMPUTER ALGORITHM 




















1. Inputs [Units] 

° Number of simulation replications, N,,,,=500. 

e Maximum simulation time, ¢,,,=150 [hour]. 

° The length of search area of X direction, / =150 
[nm]. 

° The length of search area of Y direction, I, =150 
[nm]. 

e Searcher speed, V=200 [nm/hour]. 

e Target speed, U [nm/hour]. 

e Searcher’s detection range, R=2 [nm]. 

e Searcher’s course change rate, A, [1/hour]. 

e Target’s course change rate, 4, [1/hour]. 

° The unit time of simulation, At [hour]. 

° The size of Search area, A,=150° [nm’]. 

2. Functioning of the Program 


When a new replication begins, the initial positions of 











the searcher and the target are chosen from a Uniform 





Distribution over the search area. The only difference 


between position selection logic of the searcher and the 








target is caused by the searcher’s sensor range (R). In 
order to prevent over-searching, the searcher’s initial 


position should be limited to search area (A'). 


The initial course is also drawn from a Uniform 








Distribution between 0 and 2a” . The subsequent course 


changes for searcher and target occur according to Poisson 
14 


Processes with rate 4A and 4A respectively. In particular, 





at each time step At, if Uniform_Random(0,1)</,At, then a new 
random course C,=Uniform_Random(0,1)x2z is selected for the 


searcher. The course changes for the target are determined 





in the same manner by using 4. 





When the searcher or target encounters an area boundary, 


then a random reflection occurs. After the reflection, the 





new course in radians is Uniform_Random(C,—.5, C,+.5), where C, 


is the perpendicular course from the reflection boundary. 








The parameter .5 was determined experimentally to prevent 





both “corner capture” and too many near “perpendicular 





reflections.” (see Figure 6) 
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(1) Corner effect (2) 


UniformR andor(all) UniformRandom 1 


a 
LS 


=¢, 
vant 


Perpendicular 
(4) bouncing off 
UnitormRandom- 0.5 UnitormR andom0.4 
heen aes nh, estan vine ae al Va ah | 
SAIN Need | MY RD ice =i Zl 





Figure 6. Random Movement Behavior at the Edge 
Depending on Different Parameters 


At each time step, we store the position of the 





searcher and the target, and then measure the distance 


between them. 





In order to closely approximate a continuous simulation, 


At should be small. However, too small a At means too many 





calculations, which in turn require an excessive time to 
simulate. How small a value of At is sufficient to produce 


an accurate result? 
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V=200, U=100 













I 
| 
| er 
(0). ee ge ra CE OCS EE RE SCENE DE Pee CE YOLEN Nee ea ay ag Sa eg pg cg tg eg Ne ee Z| 
I apes 
! mea a J 
ideal value 
0.8 aa 





a T 
At=0.0025 
i 











S05} -------- eee ee ee ee ee | 
x 
enc 28 eee EOE Ue Bia ee Oe 
i} i} st i 
sc a a cc rN eg ae 
| | . | | 
i i} i i 
i} i} I i 
rp ccs SL) 
if i} i} i} 
it i} i} I 
i i} i} i} 
GE CE aE PoE aE aE SE eT He ee ee a a 
i} i} i} i} 
i} i} i} i} 
i} i} i} ‘ 
if | | | 
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Figure 7. F(t) for Various Values of At 


R 
Figure 7 above shows that ee is an appropriate 


value to use. Thus the searcher moves half the distance of 





the cookie-cutter sensor’s radius at each time step. 





The author also experimented with the Poisson course 
change rate A to produce searcher and target motion most 


closely satisfying the Random Search assumptions. 


iL} 


V=200, U=0 

















Figure 8. F,(t) for Various Values of 4 


2V 
II, x, ; 


This formula implies that on the average, two course change 


As indicated in Figure 8, a good value is A= 





events occur during the time reguired for the mover to go 








from edge to edge. If A is too small, the mover has a very 


small chance to change course before bouncing off the edge. 





On the other hand, if A is too large then th mover’s 
position will be potentially limited to a small part of the 


searcher and target area (see Figure 9). 
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Figure 9. Searcher’s Movement Pattern for Various 
Values of 4 


3. Output 


Using the recommended values for the reflection 


2V 
il, x1, ’ 


the Random Search simulation produced results very close to 





parameter (0.5 radian) and course change rate (‘A= 


the expected theoretical model. 
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Figure 10. 


F,(t) for Various Target Speeds 
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In Figure 10, the pink line is plots of following 
formula: 

F,(t)=1-(1-7R?/ A,)xexp(-2RVt/A,) . (16) 

The blue line represents simulation results. The 
results indicate that this MATLAB simulation is a valid 
representation of Random Search. 

It is also clear that the movement of the target 
increases the opportunity for detection. Thus, if the 


searcher conducts a Random Search, 





target is to be stationary. 
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the best strategy for 
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IV. SIMULATION OF THE EXTENDED RANDOM SEARCH MODEL 
WHERE SEARCHER AND TARGET AREAS ARE NOT COINCIDENT 


A. DESCRIPTION OF THE GENERALIZED RANDOM SEARCH MODEL 


1. Characteristics of the Searcher and the Target 


The only difference from the previous model is that the 





search area is now not completely overlapped with the target 


area. 


B. COMPUTER ALGORITHM 




















1. Inputs [Units] 

° Number of simulation replications, N,,,,=500. 

e Maximum simulation time, ¢,,,=150 [hour]. 

° The length of search area in X direction, / =150 
[nm]. 

° The length of search area in Y direction, 1, =150 
[nm]. 

® The length of target area in X direction, /,=150 

® [nm]. 

e The: length of ‘targer area an Y directrony : 1-150 
[nm]. 

° Searcher speed, V=200 [nm/hour]. 

e Target speed, U [nm/hour]. 

e Searcher’s detection range, R=2 [nm]. 

° Searcher’s course change rate, A, [l/hour]. 

e Target’s course change rate, 4, [1/hour]. 
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R 
e The unit time of simulation, PE ag =e [hour]. 


e The size of search area, A, = 1507 [nm]. 
e The size of target area, A, =150° [nm]. 
° The size of overlap area, A,=100* [nm’*]. 
2. Functioning of the Program 


When this model was implemented in MATLAB, the faster 
mover’s At was used. For example, if the searcher’s speed 


is 200kts and the target’s speed is 100kts then, this 


Ssimulation’s At is an _ 2s hour. 


2x 200kts 





Figure 11. Target Movement Behavior for A, (V=200kts, 
U=2kts) 


The image in Figure 11 results when using the following 


parameters: ¢,, =200 hours, search area A,=150’nm’, target 
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moving area A, =100°nm’ , overlap area A, =50°nm’ , V =200kts, 








Ue Dias Re Dent, Maa 230 005 honey We ed eae 
2x 200kts © /150nmx150nm 
2x 2kts 
, and A= =0.04/hour. For both the fast searcher 





¥100nm x100nm 


and the slow target, there is an average of two course 





changes during the time the mover crosses its area. 


3. Output 





Models for the fast and slow target motion were 














introduced in Chapter . Figure 12 shows simulation results 








for several target speeds plotted against the fast and slow 


target models: 


Tes Green line : fast target assumption 

Frys: (t) = 1—exp(-2RV (A, / A, )t/ A,) (17) 
2. Pink line : slow target assumption 

Frstow(t) = (Ay / A,)(1—exp(-2RVt/ A, )) (18) 
3: Blue line : simulation results 
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V=200, U=50 V=200, U=100 V=200, U=200 























Figure 12. F(t) for Various Target Speed, Overlap Area 
A,, =100nm x 100nm 


The fast target and slow target equations can be 
combined as follows to approximate the simulation results 


for any target speeds: 


Froomvined (t) = (€XP(-EVt)) Fry (t) + = exp(-EVt)) Fras: (1) 
ue, Uw 


a l-a 


where ee [1/Vtime ]. 


JA,A.RV 


(19) 
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a decreases from 1 to 0 as the target speed U increases 
and as time goes to infinity, resulting in the fast target 


equation (17). Figure 13 shows plots of simulation results 





and Fy wneg(t) - The red line in the Figure 13 represents the 


results of equation (19). 





V=200, U=0 V=200, U=2 V=200, U=4 










































V=200, U=20 V=200, U=50 















































Figure 13. The Comparison between Simulation Results and 
F 


Tcombined (t) 


Given these results, it appears that equation (19) can 





be used as a conservative estimate (that is, a lower bound) 





of F(t) when searcher and target patrol areas are not 


identical. 
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Vv. SIMULATION OF A SYSTEMATIC RANDOM SEARCH 
AGAINST A RANDOMLY MOVING TARGET 


A. DESCRIPTION OF ‘SYSTEMATIC SEARCH’ MODEL 


Actual searchers cannot perform Exhaustive Search, but 
it is often approximated by parallel sweeps, which is 
Similar to mowing the lawn. This is called ‘Systematic 
Search’. Other forms of Systematic Search include spiral-in 


and spiral-out tracks. 














5SONM 



































(1) (2) 


Figure 14. Simple Examples of ‘Systematic Search’ 





Patterns 


Figure 14 illustrates simplified examples of a 
‘Systematic Search’ track. We examine the two types of going 
and returning ‘Systematic Search’ patterns. One is a 
‘Reverse Course Systematic Search’ in which the searcher 
goes back via the exact same track to the starting point. 
The other is a ‘Cross Course Systematic Search’ in which the 
searcher follows the track shown in Figure 14(1) and then 


goes back via the Figure 14(2) track. 
29 


B. COMPUTER ALGORITHM 




















1. Input [Units] 
° Number of simulation replications, N,,,,=1000. 
e Maximum simulation time, ¢,,,=54.72 [hour]. 
° The length of search area of X direction, 1,=148 
[nm] . 
° The length of search area of Y direction, 1, =148 
[nm] . 
° Searcher speed, V=200 [nm/hour]. 
e Target speed, U [nm/hour]. 
e Searcher’s detection range, R=2 [nm]. 
2U 
e Target course change rate, A= [1/hour] 
1x1, 
° The unit time of simulation, ae amos [hour]. 
e The searcher’s moving distance during each time 
step, iit time =VXAt=1 [nm]. 
e The size of search area, A,=148’ [nm’]. 
2. Comparison of the Searcher’s Movement Patterns 
Simulation results, illustrated in Figure 15, indicate 











very little difference between the two search paths. The 








author postulates that any reasonable Systematic Search plan 





(e.g. parallel sweeps, spiral-in, or spiral-out) will be 


effective as long as the searcher attempts 
points in the search area with the sensor on 


any point is covered twice. 
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to cover all 


time befor 
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Figure 15. Comparison Between Reverse and Cross Course 





In Figure 15, each line is defined as follows: 


Exhaustive Search (Black line) : F.(t)=(2VR1)/A, . 


Random Search (Pink line) : F.(t)=1-(—7R’/A,)xexp(—2RVt/A,) . 





Simulation Results(Blue line). 








As expected, there is very little difference between 


the reverse and the cross course’s simulation results. 
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3. Output 











V=200, U=0 V=200, U=2 V=200, U=4 
1 1 
0.5 0.5 
0 0 
Oo 10 20 30 40 50 Oo 10 20 30 40 50 50 





40 
0 


50 100 50 10 








V=200, U=8 V=200, U=12 V=200, U=14 
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Figure 16. 


F(t) for Various Target Speed 
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(2 Times 








In general, it is thought that a simulation of 





Systematic Search always performs better than Random Search. 


Figure 16 suggests that this is true. 





The measure of effectiveness used when evaluating 








search plans was “Incomplete Mean Time to Detection 


evaluated at time t ( IMTD(t) ),” defined as, 
IMTD(t ) = J d-F(x)de. IMTD(t) is used as a surrogate for 


Mean time to detection( E(T)). In fact, limIMTD(t)=E(7). We 
t—o0 





use IMTD(t) because we often cannot run the search 


simulation long enough to accurately estimate E(T) (see 





Figure 17). 


V=200, U=20 


Areal + Area2 = f..4 — F,(x))dv = E(T) = Mean time to detection 


50 
Areal = f, (1— F-(x))dy = Incomplete mean time to detection 


.. Incomplete mean time to detection < Mean tume to detection 





Time 





Figure 17. Example of ‘Incomplete Mean Time to 
Detection’ where f=50 hours. 
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We assume that the target would like to select a speed 





maximizing IMTD(t). In Figure 18, IMTD(t) is dramatically 





increased when the target speed changed from Okts to 8kts. 
And there are no significant changes between 8kts and 44kts. 
For target speeds greater than about 44kts, IMTD(t) steadily 


decreases. 


Systematic Search(V=200kts, A=148nm*148nm) 


t = 54.72 hours (2 times sweep) + 


—— f =27.36 hours (1 time sweep) | 


Incomplete Mean Time to Detection(hour) 














20 40 60 80 120 140 160 180 200 


100 
Target Speed(kts) 








Figure 18. Incomplete Mean Time to Detection for Various 
a 


Naval Postgraduate School Professor, Alan R. Washburn, 
suggested this when he wrote, “a target that wishes to avoid 
detection might actually choose to move around atU=.2V, on 


the grounds that this is enough motion to prevent an 





Exhaustive Search, but nonetheless increases the egquivalent 





searcher speed by only 1%.”4 After that speed, as Professor 


Washburn expected, the IMTD(t) decreases because the fast 





2 Alan R. Washburn, Search and detection, 4th ed. 6-3. 
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random movement of the target runs it into the searcher more 
often than away from the searcher. Figure 18 shows that even 


though the speed of target increases from Okts to 40kts, 





which is a speed enough to significantly degrade an 


Exhaustive Search, the increase of V is only 1% of searcher 


speed. 


All things considered, the author suggests that the 





optimal evasion speed for target evading a Systematic Search 


is approximately 0.05 to 0.2 of the searcher speed. 
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VI. CONCLUSIONS 


A. CONCLUSIONS AND RECOMMENDATIONS 


The main contributions of this thesis can be narrowed 
down to two results. The first result is a generalization of 


the Random Search formula which allows both searcher and 





target to move and does not require the patrol areas to be 


identical. 





The other result is that the best speed for a randomly 


moving target evading a Systematic Search ranges from 0.05 





to 0.2 of the searcher speed. 





In this thesis, the searcher always uses a cookie- 


cutter sensor. We could develop a potentially more realistic 





model by relaxing the cookie-cutter sensor assumption in 


future studies. For instance, if we assume the probability 





of target detection with sensor range R is 0.8, then 


1, if Uniform __ Random(0,1) <0.8 & range < R 


P(detection)=P, = ( otherwise 





In addition, the comparison of results when using various 
kinds of sensor types such as cookie-cutter sensor, inverse- 
cube law sensor, and triangle sensor would be worthwhile to 


study. 


Sid 
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APPENDIX 


A. MATLAB CODE 


1. Random Search Model 





Nreps=500; snumber of simulation replications 
tmax=150; Smax. Simulation time (hr) 

1x=150; Ssearch area length in x direction (nm) 
ly=150; esearch area length in y direction (nm) 
V=200; Ssearcher speed (nm/hr) 

U=50; Starget speed (nm/hr) 

R=2; Ssearcher detection range (nm) 

dt=R/ (2*max(V,U)); Sdelta t (hours) 


lams=(2*V) /sqrt (lx*ly); %searcher course change rate (1/hr) 
lamt=(2*U) /sqrt (lx*ly); S%Starget course change rate (1/hr) 



































Xs=zeros(1,tmax/dt+1); Sinitialize x-position to zero(searcher) 
Ys=Xs; Sinitialize y-position to zero(searcher) 
Cs=Xs; Sinitialize searcher course to zero 
Xt=Xs; Sinitialize x-position to zero(target) 
Yt=Xs; Sinitialize y-position to zero(target) 
Ct=Xs; Sinitialize target course to zero 
CumDet=Xs; Sinitialize cumulative detection stats 
T=0:dt:tmax; Ssimulation time vector 

A=1x*ly; search area 

for n=1:Nreps Smain simulation loop 

xs=rand* (1x-2*R) +R; Sinitial searcher and target x and y positions 








ys=rand* (ly-2*R) +R; 
xt=rand*1x; 
yt=rand*ly; 





















































cs=rand*2*pi; Sinitial searcher course 
ct=rand*2*pi; Sinitial target course 
t=0; Sset simulation time to 0 
tindex=1; Sinitiralize time index to 1 
Xs (tindex) =xs; S$save initial searcher x position 
Ys (tindex) =ys; ssave initial searcher y position 
Cs (tindex) =cs; Ssave initial searcher course 
Xt (tindex) =xt; save initial target x position 
Yt (tindex) =yt; save initial target y position 
Ct (tindex) =ct; S$save initial target course 
for t=1l:tmax/dt Sinner loop 
tindex = tindextl; Supdate simulation time index 
if rand<lams*dt; SPoisson course change rate. 
cs=rand*2*pi; SNew courses uniform (0,2*pi). 
end 
if Xs(tindex-1)<R; cs=(rand-.5);end 
if Xs(tindex-1)>(1x-R); cs=pit(rand-.5);end 
if Ys(tindex-1)<R; cs=pi/2+(rand-.5);end 
if Ys(tindex-1)>(ly-R); cs=-pi/2+(rand-.5);end 
Xs (tindex) = Xs (tindex-1)+V*dt*cos (cs); SUpdate x and y positions 
Ys(tindex) = Ys(tindex-1)+V*dt*sin(cs); 





39 


Cs (tindex) =cs; 
if rand<lamt*dt; 


ct=rand*2*pi;end 


if Xt (tindex-1)<0; ct=(rand-.5);end 

if Xt (tindex-1)>1x; ct=pit+(rand-.5);end 

if Yt (tindex-1)<0; ct=pi/2+(rand-.5);end 
if Yt (tindex-1)>ly; ct=-pi/2+(rand-.5);end 




















Xt (tindex) = Xt (tindex-1)+U*dt*cos (ct); SUpdate x and y positions 
Yt (tindex) = Yt (tindex-1)+U*dt*sin(ct); 

Ct (tindex) =ct; 

end Sinner loop (time increasing from 0 to tmax) 
CumDet = CumDet + cummax((Xs-Xt).*2 + (Ys-Yt).*2 <= R%2); 

end Souter loop (simulation replications) 


Probability=CumDet./Nreps; 
RSfun=1- (1-pi*R*R/A) *exp (-2*R*vtilde (V,U) *T/A) ; 


plot (T,CumDet/Nreps, 


"b-', T,RSfun, 'm-'), axis([0,tmax,0,1]) 
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2. ‘Patrol Areas are not Identical’ Model 

Nreps=500; 

tmax=150; 

1x=150; ssearch area length in x direction 

ly=150; sSsearch area length in y direction 

1lx_t=150; Starget area length in x direction 
ly_t=150; Starget area length in y direction 

V=200; Ssearcher speed (nm/hr) 

U=15; starget speed (nm/hr) 

R=2; Ssearcher detection range (nm) 

dt=R/ (2*max(V,U)); sdelta t (hours) 

lams=(2*V) /sqrt (1x*ly) ; Ssearcher course change rate (1/hr) 

lamt=(2*U) /sqrt (lx_t*ly_t); *searcher course change rate (1/hr) 

A=1x*ly; esearch area (As) 

Ast=100*100; Ssoverlap area 

At=1x_t*ly_t; starget area(At) 

Xs=zeros(1,tmax/dt+1); 

Ys=Xs; 

Cs=Xs; 

Xt=Xs; 

Yt=Xs; 

Ct=Xs; 

CumDet=Xs; 

T=0:dt:tmax; 

for n=1:Nreps 

xs=rand* (1x-2*R) +R; 

ys=rand* (ly-2*R) +R+50; 

xt=rand*1x_t+50; 

yt=rand*ly_t; 

cs=rand*2*pi; 

ct=rand*2*pi; 

t=0; 

tindex=1; 

Xs (tindex) =xs; Ssave initial searcher x position 
s (tindex)=ys; ssave initial searcher y position 
s(tindex)=cs; S$save initial searcher course 
t (tindex) =xt; ssave initial target x position 
t (tindex) =yt; save initial target y position 

Ct (tindex) =ct; Ssave initial target course 

for t=1l:tmax/dt 

tindex = tindextl; 

if rand<lams*dt; 

cs=rand*2*pi; 

end 

if Xs(tindex-1)<R; cs=(rand-.5);end 

if Xs(tindex-1)>(1x-R); cs=pit(rand-.5);end 

if Ys (tindex-1)<(R+50); cs=pi/2+(rand-.5);end 
if Ys (tindex-1)>(ly-R+50); cs=-pi/2+(rand-.5);end 
Xs (tindex) = Xs(tindex-1)+V*dt*cos (cs); 
Ys(tindex) = Ys(tindex-1)+V*dt*sin(cs); 

Cs (tindex) =cs; 

if rand<lamt*dt; ct=rand*2*pi;end 

if Xt (tindex-1)<50; ct=(rand-.5);end 

if Xt (tindex-1)>(1x_t+50); ct=pit+(rand-.5);end 
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if Yt (tindex-1)<0; ct=pi/2+(rand-.5);end 
if Yt (tindex-1)>ly_t; ct=-pi/2+(rand-.5);end 











Xt (tindex) = Xt (tindex-1)+U*dt*cos (ct); 
Yt (tindex) = Yt (tindex-1)+U*dt*sin(ct); 
Ct (tindex) =ct; 

end 


CumDet = CumDet + cummax((Xs-Xt).*%2 + (Ys-Yt).*2 <= R%*2); 


if n/50 == floor(n/50), n, end; 
end 
RSfun_slow=(Ast/At) * (l-exp (— (2*R*vtilde(V,U) *T) /A)); 
RSfun_fast=1-exp (- (2*R* (Ast /At) *vtilde(V,U) *T) /A); 
subplot (2,1,1) 








plot (T,CumDet/Nreps,'b-', T,RSfun_slow, 'm-',T,RSfun_fast,'g'), 


axis([0,tmax,0,1]) 

subplot (2,1,2) 

xline=[150 0 0 150 150 50 50 200 200 50 50]; 
yline=[50 50 200 200 50 50 150 150 0 O 50]; 
plot. (Xs, Ys, b="; St, Xt; re 








7XS(1),Ys(1),'yo',Xt(1),Yt(1),'yo',xline,yline,'k'),axis off 


hold on 
axis([0 1x+50 0 ly+50]), axis square 
scale=1x+t50; 














text (Xs(1)+scale/40,Ys(1),'Start', 'FontSize',10, 'FontName', 'Times"'); 
text (Xt (1)+scale/40,Yt(1),'Start', 'FontSize',10, 'FontName', 'Times"); 
hold off 
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3. Systematic Search Model 








Nreps=1000; snumber of simulation replications 
1x=148; Ssearch area length in x direction (nm) 
ly=148; esearch area length in y direction (nm) 
V=200; Ssearcher speed (nm/hr) 

U=150; Starget speed (nm/hr) 

R=2; Ssearcher detection range (nm) 

dt=R/ (2*V); $delta t (hours) 
unit_time_movement=V*dt; 

tmax=dt*10944; Smax. Simulation time (hr)=54.72 hours 


lams=(2*V) /sqrt (1x*ly 
lamt=(2*U) /sqrt (1x*ly 
L) 











hy Ssearcher course change rate (1/hr) 
); Starget course change rate (1/hr) 























Xs=zeros(1,tmax/dt+1); Sinitialize x-position to zero(searcher) 
Ys=Xs; Sinitialize y-position to zero(searcher) 
Cs=Xs; Sinitialize searcher course to zero 
Xt=Xs; Sinitialize x-position to zero(target) 
Yt=Xs; Sinitialize y-position to zero(target) 
Ct=Xs; Sinitialize target course to zero 
CumDet=Xs; Sinitialize cumulative detection stats 
CumDet_s=Xs; Sinitialize cumulative detection stats 
T=0:dt:tmax; Ssimulation time vector 
A=1x*ly; esearch area 
% define Searcher’s coordinate 
XS_S=R; Sinitial x coordinate (searcher) 
for h=1:36 
xs_column=linspace (4*h-2,4*h-2,144); xs_s=[xs_s xs_column]; 
xS_row=(4*h-1) :unit_time_movement: ((2*R) /unit_time_movement+4*h-2) ; 


XS_S=[XS_S XS_row]; 





end 
xs_f=linspace (146,146,144); 
xS_S=[xs_s xs_f]; 

for b=36:-1:1 
xs_column_b=linspace (4*b+2,4%*b+2,144); xs_s=[xs_s xs_column_b]; 
xs_row_b=(4*b+1) :-unit_time_movement: (4*b+2- 
(2*R) /unit_time_movement) ; 
XS_S=[xXS_S XS_row_b]; 











end 
xs_b=linspace (2,2,144); 
xS_S=[xs_s xs_b]; 
ys_s=2:1:146; sdefine target’s coordinate 
£67 2=12 18 
ys_row=linspace (146,146,4); ys_s=[ys_s ys_row]; 
ys_column=145:-1:2; ys_s=[ys_s ys_column]; 
ys_row=linspace(2,2,4); ys_s=[ys_s ys_row]; 
ys_column=3:1:146; ys_s=[ys_s ys_column]; 
end 
BO 2121-8 
ys_column=145:-1:2; ys_s=[ys_s ys_column]; 
ys_row=linspace(2,2,4); ys_s=[ys_s ys_row]; 
ys_column=3:1:146; ys_s=[ys_s ys_column]; 
ys_row=linspace (146,146,4); ys_s=[ys_s ys_row]; 
end 
ys_b=145:-1:2; ys_s=[ys_s ys_b]; 
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for n=1:Nreps Smain simulation loop 

xs=rand* (1x-2*R) +R; Sinitial searcher and target x and y positions 
ys=rand* (ly-2*R) +R; 

xt=rand*1x; 

yt=rand*ly; 
















































































cs=rand*2*pi; Sinitial searcher course 
ct=rand*2*pi; Sinitial target course 
t=0; Sset simulation time to 0 
tindex=1; Sinitialize time index to l 
Xs (tindex) =xs; Ssave initial searcher x position 
Ys (tindex) =ys; save initial searcher y position 
Cs (tindex) =cs; S$save initial searcher course 
Xt (tindex) =xt; save initial target x position 
Yt (tindex) =yt; Ssave initial target y position 
Ct (tindex) =ct; Ssave initial target course 
for t=1l:tmax/dt Sinner loop 

tindex = tindextl; Supdate simulation time index 

if rand<lams*dt; SPoisson course change rate. 

cs=rand*2*pi; SNew courses uniform (0,2*pi). 

end 

if Xs(tindex-1)<R; cs=(rand-.5);end 

nie a alti ndex-ty> (ix R); cs=pit+(rand-.5);end 

if Ys(tindex-1)<R; cs=pi/2+(rand-.5);end 

if Ys(tindex-1)>(ly-R); cs=-pi/2+(rand-.5);end 

Xs (tindex) = Xs(tindex-1)+V*dt*cos (cs); 

Ys(tindex) = Ys(tindex-1)+V*dt*sin(cs); 

Cs (tindex) =cs; 

if rand<lamt*dt; ct=rand*2*pi;end 

if Xt (tindex-1)<0; ct=(rand-.5);end 

if Xt (tindex-1)>1x; ct=pit+(rand-.5);end 

if Yt (tindex-1)<0; ct=pi/2+(rand-.5);end 

if Yt (tindex-1)>ly; ct=-pi/2+(rand-.5);end 

Xt (tindex) = Xt (tindex-1)+U*dt*cos (ct); 

Yt (tindex) = Yt (tindex-1)+U*dt*sin(ct); 

Ct (tindex) =ct; 
end Sinner loop (time increasing from 0 to tmax) 


CumDet = CumDet + cummax((Xs-Xt).*%2 + (Ys-Yt).*2 <= R%2); 

CumDet_s = CumDet_s + cummax((xs_s-Xt).*2 + (ys_s-Yt).*2 <= R%2); 

end Souter loop (simulation replications) 
Exhautive_Search=2*R*vtilde(V,U) *T/A; 

subplot (2,1,1) 

plot (T,CumDet/Nreps,'b-', T,CumDet_s/Nreps, 'm- 
',T,Exhautive_Search,'k'), axis([0,tmax,0,1]) 

subplot (2,1,2) 

plot (xs_s, ys_s,Xt,Yt,'r-',xs_s(1),ys_s(1),'ko',Xt(1),Yt(1),'go'),axis([0 
lx 0 ly]),axis square 
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